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Abstract 

We study equilibrium states in relativistic galactic dynamics which are de- 
scribed by stationary solutions of the Einstein- Vlasov system for collisionless mat- 
ter. We recast the equations into a regular three-dimensional system of autonomous 
first order ordinary differential equations on a bounded state space. Based on a 
dynamical systems analysis we derive new theorems that guarantee that the steady 
state solutions have finite radii and masses. 
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1 Introduction 



The Einstein- Vlasov system describes a collisionless gas of particles that only interact 
via the smoothed relativistic gravitational field they generate collectively through their 
averaged stress-energy. In this article we study the equilibrium states of the spheri- 
cally symmetric Einstein- Vlasov equations. This topic is of considerable physical and 
mathematical interest: Fackerell, Ipser, and Thorne have used these models as models 
for relativistic star clusters ^ [21 El HI El E] (see also these references for references to 
earlier literature); Martin-Garcia and Gundlach studied these models in the context of 
critical collapse [7]; Rein and Rendall approached the area from a more mathematical 
point of view jHj — this paper will be taken as the starting point for the present work. 

It is straightforward to work with particles with different masses, see e.g. pHE), but 
for simplicity we will restrict ourselves to particles with a single mass, m,Q > 0, which 
yields the same mathematical problem. A collisionless gas is characterized by a non- 
negative phase space distribution f{x^,p^) defined on the future mass shell PM of the 
tangent space TM of the spacetime, i.e., are spacetime coordinates and p^ are local 
coordinates of the four-momentum with respect to the coordinate basis x^ such that 
VyP^ = ~~ m o> P° > 0, where p^ = g illV p v ', and where g^ is the spacetime metric with 
signature (—,+,+,+); throughout this article we set c = 1 = G, where c is the speed 
of light and G is the gravitational constant. 

For static spherically symmetric models, the metric can be written as 

ds 2 = - e Mr) dt 2 + e 2A(r) dr 2 + ^2 + ^2 0#2) ( (1) 

and the Vlasov equation as 

-=L= ■ v x f - vT+^y- • v v f = o , (2) 

Vl + v 2 r 

where v is defined below. The non-negative energy-density and radial- and tangential 
pressures are given by 



P 



f(r, v)\/l + v 2 dv , (3a) 



Pr = 2 / /M)— 2 7^=^- ( 3c ) 



Prad = / f(r,v)[—) -j== , (3b) 

\x x v\ 2 dv 
r 2 y/YT 

The Vlasov equation admits a stationary solution given by f(E,L 2 ), where E and L 
are conserved quantities, interpreted as particle energy and angular momentum, per 
unit mass; they are defined by 



E := e^ r) Vl + v 2 , L 2 = r 2 v\ , (4) 

where we have followed [S] and defined v 2 as v 2 := h^u^u", where = p^/niQ and 
hfiu '■= n^riy + g^ u , where n M = g^n 11 and n M = e _At (<9j) M ; v\ is defined as v\ = v 2 — v 2 , 
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where v r = e x p r /mo and where p r is the radial coordinate component of pP. Note 
that the above objects are associated with the spatial projection of the four-velocity: 
v fj, = }i^ uU u anc i that they thus differ from the physical three- velocity, V^, which is 
defined by = r(n^ + V^) and n^V^ = 0, where V = l/\/l — V 2 ; hence, e.g., 



v 2 = V 2 /{1 - V 2 ) and E = Te^ r \ 
Using that / = f(E, L 2 ) yields: 



Prad = —^3 



00 max 

/ }{E ,L>) , E dL dE . (5a) 

r J J J \/r 2 {e- 2 ^E 2 - 1) - L 2 V ; 

/■CO /"i^ax 

/ / f(E, L 2 )y / r 2 (e~ 2 ^E 2 - 1) - L 2 dL 2 dE , (5b) 

Je^ JO 







where -L max = r 2 (e~ 2 ^E 2 — 1 



PT = ^ / f(^^)-= —= = = , (5c) 



y/r 2 (e-^EP — I) - L 2 



Let i? E (0, oo] denote the radius of support of the system, i.e., p(r) = when r > R. 
We are interested in gravitationally bound systems in equilibrium that either have finite 
radii R or possess a function p that decreases sufficiently rapidly towards infinity so 
that in both cases lim r ^#/i(r) = fj,R < oo. Because of the equilibrium assumption, it 
follows that v r of the individual particles has to satisfy lim r ^RV r = 0. We therefore 
require, since E = e^W^/l + v 2 + L 2 /r 2 , that f(E,L 2 ) = when E > E = e^, i.e., 
we assume that there exists a cut-off energy that prevents particles from escaping the 
gravitational field they create collectively. 

The basic remaining equations of the Einstein- Vlasov system are given by 

e- 2A ^-l)+l = 8vrr 2 p, (6a) 

-2A ( o^P i\ i o_„2 



\ 2 d^ ~ 1 J " 1 = ftad ' (6b) 

where £ := logr. We define the mass function, m(r), according to 

e" 2 * :=1-^, (7) 



— -=Anr 6 p. (8) 



which, together with (jba|) . yields 



Consequently, Eq. (|6b|) takes the form 



«'/' _ ,2.\ / '" , ,_ .2 



e-(- + 4vrHp rad j . (9) 

Since / is assumed to be positive when E < Eq, it follows that p > when r < R, and 
that jU is a monotonically increasing function in r (or £ = logr); it also follows that 
< m(r) < M := m(R) when < r < R, where M denotes the ADM mass. 
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In this paper we will consider distribution functions of the type 

/ = <t>(E)L 21 > (10) 

where (ft(E) is a non- negative function such that <f)(E) = when E > Eq; we also 
assume that — 1 < I < oo. 

Let us make the basic definitions 

77 := ln(£7 ) - M = - M , £:=E/E . (11) 
Rewriting in terms of ry results in 



dC V l--J ^ 7 +^«j. (12) 

Note that r? is closely related to the redshift, z, measured at the surface R: z = e v — 1. 

For distributions of the type (|1U|) the density and the radial pressure take the following 
form: 

Pr ad = ar 2l e^ l+ ^g l+3/2 (r l ), (13a) 
p = ( Z + |) ar « e (2*+4)r, [ 2ft+3/2(r? ) + e- 2 ^ i+1/2 (r7)] , (13b) 



where the constant a is given by 



9Z+3/2 3/2 + jj . 

(/ + 3/2)T(/ + 3/2) ' 



;i3c) 



here, T(-) denotes the Gamma function. The functions gi{ri) are defined by 

g t (rj):= f </>(£) [\{£ 2 - e'^f d£ . (13d) 

where, with a slight abuse of notation, <f>(E) = (f>{£); note that gi > and £ e v = 
Vl + v 2 > 1. The tangential pressure is given by pt = (I + l)p r ad> 

Throughout the paper, the ratio a = Pmd/p wm play an important role: 

a = a(r]) ■= — = ir—jj ; _ 2 < ^—57 ■ (14) 

P 3 + 2/ ^^ +3/2 + i e 2r ig l+l/2 J 3 + 2/ 

Note that this inequality is equivalent to w < 1/3, where w = p/p and where p = 
(Prad+2pT)/3 is the isotropic pressure; the limit w = 1/3, or equivalently o = 1/(3+2/), 
corresponds to the ultra-relativistic limit, i.e., radiation. 

The equations (jHJ and (|12|) form the Einstein- Vlasov system for equilibrium states; it 
is a two-dimensional non-autonomous system (recall that £ = log(r)); the system is 
closed via (|13[) which provides the relation between p rat i and p through the definition 
of gi(rj) in ()13d|) . However, note that the system Q and ()12|) is not just relevant for 
kinetic theory, but also for any fluid with an anisotropic pressure pt = (/ + l)p ra d 
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and an equation of state that relates p ra d and p (the system reduces to one for an 
isotropic perfect fluid when I = so that pt = Prad)- The fundamental assumption 
in kinetic theory that the distribution function be non-negative automatically leads 
to natural physical conditions on p, p ra d, and pt; in particular, the density and p rax i 
are positive when r < R. Since Z > — 1 it follows that p^ > as well (the limiting 
case I = — 1 describes a gas of non-colliding particles with purely radial motion, or a 
fluid with no tangential pressure). From a mathematical point of view one can regard 
the kinetic case as a special anisotropic fluid case that corresponds to a certain class 
of equations of state, implicitly determined by (j!3a|) . 1)1 3b|) and ()13d|) : however, recall 
that (j{rj) < 1/(3 + 21) in the kinetic case. 

If we assume the usual weak differentiability conditions so that T^ v = 0, where T^ u 
is the energy-momentum tensor, then we obtain the equation 

—j— = 2lp rSud + (p T3ud + p) — , (15) 

which, when (|12j) is inserted, yields the analogue of the Tolman-Oppenheimer-Volkov 
equation. From (|15|) it follows that we can write a in the form 



a 



d 



3 + 2/ + — log ^ +3/2 (r?) 



(16) 



The outline of the paper is as follows. In Sec. |2] we derive the dynamical system 
formulation: the Einstein- Vlasov system is reformulated as an autonomous system of 
differential equations on a bounded state space. This dynamical system is subsequently 
analyzed in Sec. 01 The results are then interpreted in the standard physical variables 
to yield our main theorems in Sec. 0J these theorems formulate various criteria that 
guarantee solutions of the Einstein- Vlasov system to have finite masses and radii; hereby 
we generalize earlier work (Hj. We conclude with a brief discussion in Sec. 03 In 
Appendix E] we investigate various differentiability conditions that are needed in the 
paper and derive properties of the prominent function cr(r]). In Appendix[Blwe discuss 
existence and uniqueness of solutions. In Appendix O we perform a center manifold 
analysis, which allows us to sharpen one of the theorems when certain asymptotic 
differentiability conditions are imposed. Finally, in Appendix D we present an extension 
of the monotonicity principle given in jlUj . 



2 Dynamical systems formulation 

In the following we reformulate the Einstein- Vlasov system as an autonomous dynam- 
ical system on a bounded state space. As a first step we introduce dimensionless 
asymptotic homology invariant variables jllj 

47rr P T /-in 
u= , v = — ^ — -: (17) 

m a (i-2my v ; 

the variable v is not to be confused with the velocity variable of Section ^ In the 
variables (n, v, rf) the Einstein- Vlasov system as given by Q and (|12|) takes the form 



2 DYNAMICAL SYSTEMS FORMULATION 



6 



of an autonomous system, 

du 

dv 

di 

drj 

where 



[3 + 21 -u - (l + a — a ') h] (18a) 
[{u- l)(l + 2av) + a'h] (18b) 
-ah, (18c) 



h := (1 + au)v ; (18d) 
note that henceforth a' denotes da/drj. 

In many cases, depending on the choice of the distribution function (f)(£) , cr'(r]) is a 
rather complicated function; it is thus convenient to replace 77 by a variable u; that 
is adapted to the choice of (ft, so that a' as a function of cj has a simpler form. In 
addition, the degree of differentiability of a' at 77 = can be improved; for instance, the 
function <r' = C\ + C^y^ is not C 1 at 77 = 0; however, by choosing w = y/77 we obtain 
a' = C\ + C2W, which is smooth. 

We introduce the variable to by choosing oj{rj) according to the following rules: uj(rj) 
is smooth on (0, 00) and strictly monotonically increasing, i.e., du/dr] > 0; oj{rj) > 
when rj > 0, and uj — > when 77 — ^ 0. Replacing r/ by u; in the system (|18|) leads to 

duJ dto , „ , , „ dlog a; 

— = -a — h = -uFh where F := a — ^— ; JTO ) 
ct^ otr/ «77 

we assume that i 7 = a dlogco/dr] is bounded for all r? and in the limit rj — > 0. The 
specific choice of a; depends on the type of distribution function one is considering: as 
we will see below, besides the trivial u = 77, one useful choice is 

uj = k (e br i - l) , (19) 

where k > and b > are constants. In this case, doj/dij = kbe bri so that duj/dr] — > const 
as rj — > 0; moreover, F = ba(l + k/to), which is bounded in the limit 77 — ^ when a' 
is bounded. In general, it turns out that it is often advantageous to choose u to be of 
the form 

uj = k(ip( V ) - 1) , (20) 
where /c is a positive constant and 93(77) satisfies cp > 1, dip/dr} > 0, and </?(0) = 1. 
The last step is to introduce bounded variables (U,V,Q) based on (u,v,uj): 

u --=tt- v -=rir- n ^v^-- < 21 > 

1 + U I + V 1 + OJ 

furthermore, let us introduce a new independent variable A according to 

^ = (l-U)- 1 (l-V)- 1 . (22) 
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This leads to the following dynamical system on the bounded state space of the variables 

(u,v,ny. 

^L = U(l- U) [[3 + 21 - (4 + 2l)U]{l - V) - (1 + a - a') H] , (23a) 
CLa 

= V(l - V) [(2U — 1)(1 — V + 2aV) +a'H], (23b) 

CLa 

^ = -0(1-0) FH, (23c) 
dX 

where 

H := (1 - U + aU)V ; (23d) 

the quantities a, da/dr], F are now to be regarded as functions of O. 

The dynamical system (|23|) is polynomial in U and V, so that the system extends 
to the boundaries U = 0, U = 1, V = 0, V = 1. The state space is thus given 
by (U, V, O) G [0, l] 2 x (0,1). The limit a — > corresponds to the Newtonian limit, 
since the system (1231) then reduces to the system of equations for the Newtonian case; 
the Newtonian equations (for the polytropes) are also recovered for O — > (which 
is connected with a — > 0, see ()24|) below); however, in order for that limit to exist, 
lim^o a ' must be defined. Hence, before we can proceed to analyze the global dynamics 
of the system (|2*3*)l we must discuss the properties of a' , in particular its degree of 
differentiability and its limiting behavior as ij — > 0. 

Throughout the paper we will assume that cr'(r]) is a continuous function on r] G (0, oo), 
or, equivalently, that o' as a function of O is in C°(0, 1). In Appendix we discuss 
in detail which assumptions on the distribution function are required to achieve a 
continuous function a 1 , here we merely state the basic facts: we consider a distribution 
function such that 4>{£) satisfies: 4>{£) = for £ > 1; (f> € C°(0, 1), and 0(f) < 
const for some e > in a neighborhood of £ = 1. If — 1 < Z < — 1/2 we assume, 

in addition, that is Holder continuous with index a > —I — 1/2. 

By assumption, a' is continuous for rj > 0. lim^o c'^) need not exist; however, from 
now on we will assume that a' has upper and lower bounds in the limit r] — > 0, i.e., 
fci < a' < &2 and hence k\r] < a < kif] as r\ — > (where fei, &2 are positive constants). 
A direct consequence is 

lim a(rj) = . (24) 

Higher differentiability of a' requires higher differentiability of <j){£). 

In Appendix^we discuss the 'asymptotic polytropic case'. For asymptotic polytropes 
we find in addition to (l24l that 



a' := lim a' = , (25) 

u r?->o 1 + Z + n v ; 

where the constants I and n satisfy — 1 < I < oo and 1/2 < n < oo (the latter condition 
is motivated from kinetic theory: n > 1/2 implies n — 3/2 > —1, cf. (|26|) . so that 
5j(r/) exits). The bounds on Z and n imply that < a' < 2. A statement equivalent 
to (|25ft is that a' is in C°[0, oo) in the asymptotically polytropic case. As we discuss in 
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Appendix El the asymptotically polytropic case requires a distribution function <j){£) 
with asymptotically polytropic behavior, i.e., 

</>(£) = <M1 " £) n ~ 3/2 + (26) 
as £ — > 1 (with n > 1/2); 0_ > is a constant. 

Let us conclude this section with an example: the relativistic polytropes. This class 
of isotropic pressure models (i.e., / = 0) is characterized by the following equation of 
state: 

p = Kp 1 ^ +p/{T -1) cp = [x 1/r + x]/(T - 1) , (27) 

where T, K are constants and c = [(r — 1)-K~]~ r ' ( r ~ 1 ) , x = cp. In it was shown that 

u = x (r-i)/r = i {e v_ lh (28) 

which is of the form to = k(e br) — 1) with k = 1/F and 6 = 1, leads to 

* = Cr-i)n, % = = ^ ((r ~ 1)0 + ~ n) ' (29) 

and a dynamical system with purely polynomial right hand sides. 

In [U [3] the kinematic interpretation of the following class of equations of state was 
developed: p = Kp 1 ^ 1 +p/(Y2 — 1) (also denoted as relativistic polytropes in the 
given reference); when I\ = I^ = T we recover the previous example. The associ- 
ated distribution functions are given in [1], but reasonably simple expressions are only 
obtained for r"i = T2 = T and for special values of T. For instance, T = 5/4 yields 
4> = <j)_S-S(l-£ 2 fl 2 . 



3 Dynamical systems analysis 

In Appendix [B] we establish (global) existence and uniqueness of solutions of the dy- 
namical system (|2*S|) on the state space X = (U, V, Q) £ [0, l] 2 x (0, 1), where we assume 
that a' be continuous on (0, 1). In this section, we analyze the global dynamics of the 
flow of the dynamical system on X. 

The key ingredient in the dynamical systems analysis is the fact that f2 is strictly 
decreasing in the interior of the state space X, which entails that the a- and u-limit sets 
of orbits (e.g., fixed points) reside on the boundaries of X. Standard local dynamical 
systems analysis requires the dynamical system to be at least C 1 — in that case the 
system can be linearized at the fixed points. The dynamical system (|23|) is C 1 on 
the state space [0, l] 2 x (0,1), if a' is in C 1 (0, 1). However, we will not make this 
smoothness assumption on a'; for the following we merely assume a' to be in C°(0, 1); 
in our proofs, monotonicity arguments will replace the standard arguments from local 
dynamical systems analysis. 

The inclusion of the boundary 17 = in the dynamical systems analysis depends on 
the differentiability properties of a' at f2 = 0. If a 1 is in C°[0, 1) (instead of merely in 
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C°(0, 1)), then the dynamical system continuously extends to the state space [0, l] 2 x 
[0, 1); this is the case for asymptotically polytropic distributions. If a' (as a function 
of $7) is in C^O, 1), then the dynamical system smoothly extends to [0, l] 2 x [0,1). 
Although we will not make these assumptions in the following, it is useful to consider 
the case of smooth a' to obtain a more intuitive feeling for the properties of the flow of 
the dynamical system. Figure ^ depicts the state space X and its extension to £1 = 0; 
for figures describing the dynamics on f2 = when I = 0, see the appendix in |12j . 
Tabled hsts the fixed points and the eigenvalues of the linearizations of the dynamical 
system. 




Figure 1: The state space and fixed points; note that -B4 only exists for certain a' (or 
equivalently n) values and that an additional line of fixed points appears when a' Q = 1, 
see Tabled 



Fix points 


U 


V 


n 


Eigenvalues 


Restrictions 


Li 
L 2 
L3 
U 
Bi 

B 2 
B3 
B A 


1 

3+22 
4+22 



U 



1 

3+l-n 






1 

1 

1 
1 

2(1+0 


n 
n 
n 








1; 1; 

3+22. 1+2. n 
4+22' 2+2' u 

-1; 3 + 21; 

-(l-E/b); A 3 ; 

2(2+2)^,-1. a' , - Fn 
l+a' ' HV > 2+a' 

-(1-a'o); -a' ; -F 
0; 0; 
a[l±(5 + 3/-n) v / 6];A 3 


a' = l 
3 + /-n>0 


2(2-n) 


3+22-2(7,', 



Table 1: Fixed points and associated eigenvalues; the third eigenvalue A3 always has 
the same form for all fixed points on f2 = 0, namely A3 = — -Fo^(l — U), Fq = F(0); 

« = 2 % + {n+i- % ' b = 1 + 221 + 3312 + 8/3 + 2 ( U + l3l ) n - ( 7 + 8 0« 2 > wh ere U and V 
take the relevant fixed point values. 



In the following we investigate the global dynamics of the dynamical system on the state 
space [0, l] 2 x (0,1). The most useful tool for our analysis will be the monotonicity 
principle, see Appendix El 
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Lemma 3.1. The uj -limit of every interior orbit is located on £1 = 0. 



Proof. From equation (|23c[) we see that fHs a strictly monotonically decreasing function 
on the state space [0, l] 2 x (0, 1) except for on the side face V = 0, where f2 = const. It 
follows from the monotonicity principle that the w-limit of every interior orbit is located 
either on O, = or V = 0. We prove the lemma by excluding the latter possibility: 

Suppose that there exists an orbit 7 that possesses an w-limit point on V = (with 
U 7^ and U 7^ 1). Then a point P on L2 must be in a; (7) as well, as follows from 
considering the structure of the flow on the side face V = 0. Since dV/dX > in a 
neighborhood U of L2, the orbit 7 cannot remain in U, but must leave and re-enter U 
infinitely many times. The boundedness of the r.h.s. of the dynamical system implies 
that the return time is bounded below and that O must decrease by at least AQ > 
during that time. However, this is in contradiction with the fact that P is an cj-limit 
of 7. Using similar (simpler) arguments we can also exclude that L\ and L3 are eligible 
cj-limits. Therefore, the side face V = does not contain w-limit points of any interior 
orbit, which leaves O = as the only attr actor. □ 

Remark. The proof of Lemma 13.11 does not require the assumption that the r.h.s. of 
the dynamical system be C 1 on [0, l] 2 x (0, 1) — we have merely assumed continuity 
of a' on (0, 1). However, if C 1 is assumed, the line of argument can be simplified by 
noting that L\, L2, L3 consist of transversally hyperbolic fixed points: points on L\ 
are sources, points on L2 and L3 are saddles. 

Lemma 3.2. The a-limit set of an interior orbit is either located on L\, L2, or Vt = 1. 
A two -parameter set of orbits originates from L\, a one-parameter set from L2. 



Proof. The monotonicity of £1 implies that the a-limit of an interior orbit must lie either 
on Q = 1 or V = 0. Points on L3 cannot acts as a-limits: the line of argument resembles 
the proof of Lemma l3.11 where we use that dU/dX > and dV/dX < in a neighborhood 
of L3. In a neighborhood of L\ we observe that dU/dX < and dV/dX > 0; hence, L\ 
is a source for a two-parameter family of orbits. Orbits originating from L2 correspond 
to regular solutions of the Einstein- Vlasov equations, as we will see in Section @J The 
existence of a one-parameter set of such solutions has been proved in ^5]. D 

Proposition 3.3. Suppose there exists 77 > 0, arbitrarily small, such that o~(rj) satisfies 

17 + 4/ 

a' > + cr 30) 

-4 + 2/4 + 2/ v ; 

for all 7] < i). Then the uj -limit of every interior orbit lies on (Q = 0) H (V = 1). 

Proof. Let be the value of O associated with fj. As the first step in the proof we 
show that the function Z defined by 

U , V 



l-UW-V 



is strictly monotonically increasing on (0, l) 2 x (0, n) if and only if condition (jSOJ) holds. 



3 DYNAMICAL SYSTEMS ANALYSIS 11 



A straightforward computation leads to 

Z~ x — = [2(2 + l)a - 1 - (7 + 4J)ff][V(l -U) + aUV] + A > A > ; (32) 

here, we have used that A := J7[2(l + l)(l - V) + 2(3 + 2l)Va(l + a)} > and the 
assumption (J3UJ). Conversely, if cLZ/cZA > for all (f7, V) G (0, l) 2 (and < A), 
then dSOJ) follows by letting r/ -» in (l32l). 



On the side faces {(E/, V) | f7 = 0, V < 1} and {(L7, V) | V = 0, U < 1} the function Z 
assumes its infimum: Z = 0. The monotonicity principle, cf. Appendix [DJ states that 
w-limits cannot be located on this set. As a consequence, the w-limit of every orbit 
of the considered set (0, l) 2 x (0,O) must be contained on (U = 1) or (V = 1). Since 
— > for all interior orbits, we find that the w-limit of every orbit in (0, l) 3 lies on 
(Q = 0) Pi (U = 1) or (fi = 0) Pi (V = 1) (or the union). It remains to prove that the 
set (f2 = 0) D (V = 1) is the exclusive attractor. 

Consider Eq. (l23bl) : dV/dX oc [{2U - 1)(1 - V + 2aV) + a'H]. In the case U > 1/2, 
dV/dX is manifestly positive. In the case U < 1/2, positivity is guaranteed under the 
condition 

1-2U l-V + 2aV 
> l-[/ + aC7 V ' (33) 

Since (1 - 2f/)/(l - U + <rC/) < 1 for all U G [0, 1/2), a sufficient condition for to 
hold is a' > v^ 1 + 2c, which, by assumption (|3U|) . is satisfied if i? > 2(2 + /)(1 — <r) _1 . 
Therefore, there exists an e > such that dV/dX is positive on S := [0, 1] x (1 — e, 1) x 
(0,e). Consequently, S is a future invariant set of the dynamical system and V is a 
monotonic function on S; by the monotonicity principle, the a>limit of every orbit in 
S must lie on (Q = 0) n (V = 1). 

Now suppose that there exists an orbit 7 whose w-limit a; (7) contains a point on the 
set (p, = 0) fl (U = 1) (with 7/1). The flow on this set is known explicitly; it follows 
that the point (1, 1,0) is in a; (7). Consequently, for all e > there exists X £ so that 
V(A e ) > 1 — e and 0(A e ) < e. For sufficiently small e we find 7(A e ) G 5, whereby 
7(A) G 5 VA > A e , since 5 is future invariant, and consequently V(X) — > 1 as A — > cxd. 
This is a contradiction to the assumption about the cj-limit of 7. Accordingly, the 
w-limit of every orbit must he on (D, = 0) fl (V = 1). □ 

We now consider the surface 

P(u,v) :=2u + v- 2(3 + 2/) =0 (34) 

in the state space of the variables (u,v,uj); equivalently, P = generates the surface 
[4(2 + Z) - (9 + 4l)V]U + (7 + Al)V - 2(3 + 2/) = in the state space X = {(17, V, Q) G 
[0, l] 2 x (0, 1)}. However, since u and v are finite on P = and the associated equations 
are simpler, it is preferable to perform the following analysis in the variables (u, v,u). 

Lemma 3.4. Suppose that there exists f] > 0, such that o~(j]) satisfies 

a > , 1 - + k a (35) 
~ 2(3 + 21) y 1 
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for all rj <fj, where k = 3 + 21 + ~ if I > — §, and k = 2 if — 1 < I < — |. T/ien 



p=o 



> 



(36) 



wz/ien u> < w fit 7^ 0, t> 7^ Oj; Ziere cD is the value of uj associated with fj. 
Proof. Using the system (|18p we obtain by a straightforward calculation 

- 1 + 2sa' - 2<jV + a (2(1 - 3u)u + 2s(-2 + (2 + a')u)) 



dP_ 



P=0 



where s = 3 + 21. It follows that 



dP 



> 44> a' > h k(u,ri)a , 

p=o 2s 



where 



-As{u - 1) + it(-3 + 2(3 + a)u) 



2s(l + au) 

Assume that I > — | , so that s > | , and define k = s + \ . We have 



(37) 

(38) 
(39) 



2s(l + au) k(u,n) - k = -(s - u)[2s - 3 + 2u(3 + a)] - s(2s - l)au < , (40) 

where we have used that < u < s on P = 0. Accordingly, k > k(u,r]) on P = 0; 
therefore, from (|35|) . we conclude that (|3"5|) is a sufficient condition for (|3l))) to hold. 
Incidentally, note that in (|40j) equality holds in the limit n — > (a — > 0) and u — > s; 
hence, k = s + ^ is the smallest constant such that (j35j) implies (|36|) . The proof in the 
case I < — j is analogous. □ 



Proposition 3.5. Suppose that there exists fj > 0, such that a{rf) satisfies 



a' > 



2(3 + 2/) 



/c a 



(41) 



/or a// J] <fj, where /c = 3 + 2/ + ^ if I ~> —\, and k = 2 if — 1 < I < — |. Consider an 
orbit 7 whose a-limit is a point on L\ or L2 with £1 < Cl (ij < fj). Then the uo-limit of 
7 is located on (fl = 0) H (V = 1). 

Proof. Consider the region P of the state space defined by P > and Q < Cl; this set is 
future invariant by Lemma 13.41 On the set V, V is a strictly monotonically increasing 
function. To see this, recall that dV/dX is manifestly positive when U > 1/2, and from 
manipulating (|33|) that 

1-u 



a' > 



-y- 1 + 2a 



(42) 



1 + au 

is a sufficient condition for dV/dX > on U G [0,1/2) (V 7^ 1). By assumption 
we have a' > 2 (3+?i) + ka > 2(3+21) + > therefore, condition (|l2*|) is satisfied if 
(1 — -u) -1 (l + o~u)v > 2(3 + 2/). It is straightforward to show that this relation holds 
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for (u,v) such that P(u,v) = 0, from which follows that it is satisfied for (u,v) such 
that P(u, v) > as well. This establishes the claim that V is monotonic on V. 

Since V is future invariant and V is a monotonically increasing function on V, the 
monotonicity principle in conjunction with Lemma 13 . 1 1 implies that the w-limit of every 
orbit in V resides on (Q = 0) n (V = 1). Orbits that originate from L\ satisfy P > 
initially, i.e., these orbits are contained in V . For an orbit 7 whose a-limit is a point 
on L2 we have P — > as A - > — 00. However, a simple computation shows that the 
unstable manifold of every point on L2 (with VL < ti) lies in V . The orbit 7 thus lies 
in V for all A and must possess an w-limit on (fi = 0) fl (V = 1). This concludes the 
proof of the proposition. □ 

Remark. In the asymptotically polytropic case, ctq = lim^o a' exists: o~' = (l+l+n)^ 1 . 
Therefore, in the asymptotically polytropic case, the assumptions of Proposition 13.31 
see (|3U[) . are satisfied if 

n<3 + l. (jSOI) 

Furthermore, 

n < 5 + 31 gU) 

is sufficient for condition 1)41(1 in Proposition 13.51 to hold: namely, if n < 5 + 3/ (a' > 
(6 + 4/) _1 ), then there exists fj such that a' > (6 + 4/) _1 + fccr for all rj < fj. 

In the asymptotically polytropic case, cr'(r/) is continuous on 77 S [0, 00). If a' (as a 
function of fi) is even C 1 [0, 1), then the dynamical system extends smoothly to the state 
space [0, l] 2 x [0, 1). Smoothness of o-'(ry) requires smoothness assumptions on (f>{£), 
which we refrain from discussing here; let us merely note that, e.g., polytropic functions 
4>(£) = 4>-(l — £) n ~ 3 / 2 , and linear combinations thereof, satisfies the requirements. If 
(j'{ri) is C 1 (or, more generally, if there exists a variable ui such that a' is C 1 as a function 
of u), then Proposition 13.31 can be improved: 

Proposition 3.6. Assume that there exists a variable ui such that ct'^aa is C 1 ^, e) 
(for some e > 0). If 

Op > ^_^_2i ' or ' e Q u i va l en tly> n < 3 + / , (43) 

then the id-limit of every interior orbit lies on (O = 0) D (V = 1). 

Remark. Recall that Q, can be chosen to be a smooth function of u; hence a' is C 1 as a 
function of id if and only if it is C 1 as a function of Q. In the simplest case, a' is C 1 ^, e) 
as a function of rj and thus automatically satisfies the assumptions of the proposition. 

Proof. The proof of the proposition is based on center manifold theory; it is given in 
Appendix O □ 



4 Physical interpretation 

We define a regular solution of the static Einstein- Vlasov equations as a solution that 
possesses a regular 'potential' fj,. As r — ► 0, fj,{r) converges to a central value /i c ; 
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equivalently, rj(r) —* ry c as r — > 0. This implies that the limits lim r _^o( r_2 V) an d 
lim r ^o( r ~ 2 Vrad) exist, as follows from ((T3|) , Accordingly, p — ► oo as r — ► 0, if / < 0; 
p — > const, if / = 0; p — ► 0, if I > 0. Hence, regular solutions are characterized by a 
regular metric; however, p and p ra( j are only regular when / > 0. 

For regular solutions, the variables U and V satisfy U — ► (3 + 2/)/(4 + 2/) and V — ► 
as r — > 0; this is a direct consequence of (fTTj) (where we use dm/dr = 4itr 2 p and the 
behavior of p as r — > 0). In the dynamical systems formulation, regular solutions thus 
correspond to orbits whose a-limit is a point on the line Li- 
Remark. Orbits whose a-limit is a point on L\ correspond to solutions with a negative 
mass singularity at the center. For a discussion of such solutions (in the Newtonian 
case) see [T2*] . 

Lemma 4.1. An orbit in the state space with co-limit on (Q, = 0)n(V = 1) corresponds 
to a solution of the Einstein- Vlasov equations with finite radius and mass. 

Proof. Integration of the defining equation for A, dlogr/dX = (1 — J7)(l — V), leads to 



r = tq exp 



A 



[l-U(X)){l-V(X))dX 



(44) 



where Ao and ro are such that (U,V,Q)(Xq) corresponds to (m(ro),rj(ro),ro). In order 
to prove that the solution possesses a finite radius we show that r^i?<ooasA^oo. 



Define 5U = 1 — U and 5V = 1 — V. Using these variables, Eq. (|23b|) takes the form 

— — = (5V)V [5U(a'(l - a)V - 25V - 4aV) + 5V + aV(2 + a')) . (45) 
dX 

By assumption, a — > and V — > 1 (5V — > 0) as A — * oo; therefore, for sufficiently large 
A, 

— 77- > (5C/<5F<t'(1 - e) (46) 
dX 

where e > is small. Inserting (|46[) into (|44|l leads to 

/ 5U5VdX= 5USV(-^-) d{5V)< ,-d(5V). (47) 

Jx Jsv(X) V dX ) Jsv(X) (1 - e) <r 



This integral remains finite in the limit A — > oo (<5V — > 0), since cr' > const > for 
small 7] (and thus for small SV); hence r — > i? < oo in the limit A — > oo. Finiteness of 
the mass follows from the finiteness of the radius. □ 



In combination with Lemma [4. 11 the results of Section |5] imply theorems that formulate 
criteria guaranteeing finiteness of solutions of the Einstein- Vlasov equations. 

Let us first recapitulate the basic assumptions for the theorems: we consider a distri- 
bution function of the type / = <p(E)L 21 with I > —1. Assume that there exists Eq > 
such that 4>{E) = for E > E , (j) G C°(0,E o ), and 4>(E) < const (E - E)~ 1+e for 
some e > in a neighborhood of E = Eq. If — 1 < I < — 1/2 we assume, in addition, 
that (j) is Holder continuous with index a > —I — 1/2. 
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Theorem 4.2. Suppose there exists fj > 0, arbitrarily small, such that cr(rj) satisfies 

17 + 4/ 

" a 4T2i + 4T2i" (48) 

for all rj < fj. Then every solution of the static Einstein- Vlasov equations has finite 
radius and mass. 

Remark. In the asymptotically polytropic case, provided that cr'\ v (u) £ C 1 [0, 1), the 
condition a' > 1/(4 + 2/) n < 3 + /) is already sufficient. 

Remark. Theorem 14.21 is a natural generalization of jjSl Theorem 3.1], where it is as- 
sumed, in our notation, that / > —1/2, n + I > 1, and n + / < 3. 

Proof. The theorem follows immediately from Proposition 13.31 and Lemma 14.11 □ 
Theorem 4.3. Suppose that there exists fj > 0, such that a(rj) satisfies 

for all rf < fj, where k = 3 + 21 + ^ if I > — |, and k = 2 if — 1 < I < — |. T/ien 
every regular solution of the static Einstein- Vlasov equation whose central value r/ c of 
the potential is smaller than fj, i.e., r\ c < fj, has finite radius and mass. 

Remark. Finiteness of the radius and the mass also holds for solutions whose a-limit 
lies on L\, which correspond to solutions that possess a negative mass singularity. 

Proof. The theorem follows by combining Proposition 13.51 and Lemma 14.11 □ 

Remark. In the asymptotically polytropic case, the condition n < 5 + 31 a' Q > 
1/(6 + 4/)) is sufficient for (|49|) to hold for some fj, see CT . Therefore, if n < 5 + 3/, 
then all regular solutions with sufficiently small ij c have finite radii. 



5 Concluding remarks 

In this paper we have recast the stationary spherically symmetric Einstein- Vlasov equa- 
tions for distribution functions of the form / = cj)(E)L 21 into a three-dimensional dy- 
namical system with compact closure. This has allowed us to derive new theorems 
formulating conditions under which solutions have finite radii and masses; our the- 
orems naturally extend previous work, see e.g., jSj El and the Newtonian analysis 
in [Hj. 

In the present paper we have focused on distributions that can be bounded by asymp- 
totic polytropes in the low r] regime. Additional results can be obtained if one also 
makes restrictions that lead to controlled behavior in the large r\ regime; e.g., one can 
investigate classes of distributions that yield asymptotically linear relations between 
the pressures and p, which in turn may make it possible to obtain theorems about 
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mass-radius diagrams of the type given in for perfect fluids with asymptotically 
linear equations of state in the high pressure regime. 

The theorems presented in this article may be less sharp than one would wish, since 
the natural kinetic assumptions imply non-trivial restrictions on the effective implicit 
equation of state that have not been investigated here. The connection between the 
distribution function and the effective equation of state is given via the integrals of 
Section ^ it is an interesting question whether it is possible obtain more comprehen- 
sive information about this relation, e.g., by finding classes of functions <f>(E) that 
lead to certain types of equations of state. Simple results in this direction are due to 
Fackerell [I]: it could be excluded that models with constant energy density possess a 
kinetic interpretation. In this context one should note that the assumption I 7^ im- 
plies that the density and the pressures depend on both r and n; however, the function 
a appearing on the r.h.s. of the equations is a function of rj alone. Hence the mathemat- 
ical correspondence between fluids with equations of state and kinetic models amounts 
to a correspondence in o~(r]) alone — which equations of state yield c{rj) with kinetic 
interpretations? 



A Differentiability and asymptotic behavior 

The r.h. sides of the dynamical systems ()18|) and (|23|) contain the function o-{rf) and its 
derivative. To analyze the differentiability properties of these function we investigate 
the key quantities 

gM= T m\W 2 - e ~ 2V )]' > ( 50 ) 



on which a(n) is built, see ()13d|) and (fTl)) . 

Lemma A.l. Let <f> : (0, 1) 3 £ 1— ► R be non-negative. Assume that cf) £ C k (0, 1) and 
that |(1 - £) fe (fe) (£)| < const (1 - £)~ 1+e for some e > in a neighborhood of £ = 1. 
Then, for i > 0, g^rj) £ C k+1 (0, 00) and 

^ 9i (v) = ie-^ 9i -x (i>0), ±g ( v ) = e-^(e^). (51) 
an drj 

In the case —1 < i < 0, giirj) £ C fc (0,oo); however, if the additional Holder continuity 
condition (1 - £) k ^ k \£) £ C a (0, 1) with a > —i holds, then again giirj) £ C fc+1 (0,oo). 

Proof. The idea of the proof is to transform the function gi to a form that corresponds 
to Equation (Al) of ^1] and to use the results of 14, Appendix A]. Employing an 
adapted variable u> instead of 77, cf. (|20|). the function gi can be regarded as a function 

gi(u>); we set 

u = l-e- r >; (52) 

proving differentiability in 77 amounts to proving differentiability in uj. Note that 
< to < 1 since < rj < 00. Replacing £ by x = o;~ 1 (l — £) in the integral then 
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transforms (jHOjl to 



gi(u) = / c/)(l-ujx)G l (u,UJx)(l-x) l dx, (53) 

JO 

where G{u,r) := 1 — | — ^ is a smooth, positive, and bounded function on (u, r) £ 
[0, 1) x [0, 1); in particular, x i— > G*(u;, ux) is smooth, positive, and bounded on x £ [0, 1] 
for all uj. Define 

<ji(v,u) :=uj i+l [ <P{l-ujx)G i (u,ujx)(l-x) i dx; (54) 



J o 

evidently, <7i(u;) = gi(u,uj). Introduce fa according to fa{v,r) = <p(l — t)G 1 (u,t) . By 
assumption, r i— ► fa(u,r) is measurable, non-negative, bounded on compact subsets of 
(0,1), and C k ; further more, for small r, \r k d k fa(u,T)\ < const r 1+e for some e > 0. 
By using fa, Eq. (jHljl becomes 

gi(v, uj) = uj 1+1 / fa{y , ux) (1 — x) 1 dx . (55) 
J o 

Therefore, w i— > gn(u,uj) is of the same form as Equation (Al) of ^3]. Using the results 
of Appendix A] we conclude that uj h-> g^v^u) is C k+1 (where the additional 
requirement that the function (1 — £) k (f>^ be Holder continuous enters in the case 
— 1 < i < 0). Since v i— > gi(y,uj) is smooth we find that w i— > <7j(cj,u;) = <7j(o>) is C fc+1 , 
as claimed. A straightforward computation based on the relation 

-^9i(v) = (d u gi)(uj,uj) + (d u gi)(u,u) (56) 

yields En. (|5T)t. □ 

Remark. In the case — 1 < i < the derivative of is given by 

^(r?) = w i e- (m) V(e"' ? ) + ie" 2 V / #(r/, x)(l - x^dx, (57) 

Jo 

where 5 = a + i > and u = 1 — e _?? ; since is Holder continuous with index a, the 
function *(ry, x) = ^(e" , » + w(l - a;)) (e - " + f (1 - x)) i_1 - ^(e^Ke^)*- 1 (l-x)" a 
is bounded and the integral is well-defined. The derivation of this formula is analogous 
to the considerations of j!41 Appendix A]; we omit the details here. 

Lemma lA. II implies the following corollary: 

Corollary A. 2. Let <f> : (0, 1) 3 S i— ► K be non-negative. Assume that <f> £ C k (0, 1) and 
that |(1 - £) k ^ k \£)\ < const (1 - £)~ 1+t for some e > in a neighborhood of £ = 1. 
Then 

cN +fc+1 (0,oo), (58) 

where [i] denotes the largest integer less or equal to i. For i ^ N, i/ in addition (1 — 
£)V (fe) (£) e C Q (0,1) with a > 1 + [*] -i, taen ^(77) e CW +/c+2 (0, 00). 
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Proof. Suppose that N ^ i > 0. Applying (|51j) iteratively [i] times we find that 
dfigi/drp} oc g^uy Since > 0, <?j_[j] is in C k+1 by Lemma [A. It accordingly, 

gi G C^ +k+1 . The case i S N is analogous: i = [i], d % gijdrf oc go, and 50 £ C fc+1 . If in 
addition (1 - £) k ^ k \£) e C a (0, 1) with a > 1 + [t] - i (i N), then ( ^ +1 5i /d7 ? [ i ]+ 1 oc 
ffi— W— 1) which is C fc+1 by assumption; from this the claim follows. □ 



Recall from (|16|) that 
Accordingly, 



3 + 21 + -^- log g l+3/2 (v) 



(59) 



do- o ci 2 



2 1 °Sft+3/2('7) • (60) 
dr) dr] z ' 

Under the assumptions of Corollary IA.2I we see that g l+3/2 is in C^ l+3 ^ +k+1 so that 
da/dr] is in cl>+3/2]+k-l ( an d in c['+ 3 / 2 ]+ fc when the Holder condition is satisfied). In 
particular, if I > —1/2, then the assumptions <fi G C°(0, 1), |0(f )| < const (1 — £)~ 1+e 
for some e > as £ — ► 1, i.e., = in Corollary IA.21 suffice to obtain gi +3 /2 S C 2 and 
da/dr] £ C°; in the case Z < —1/2, the additional Holder condition must be imposed. 
Higher differentiability of (f> (larger k) leads to higher differentiability of da/drj. 

Definition. A function gi(rj), i > — 1, is said to possess asymptotically polytropic 
behavior if there exists R B n > 1/2 such that 

gi {ri) = const rf +n ~ 1 l 2 (1 + o(l)) (61) 

as T) — > 0. VFe require the derivatives of gi, if they exist, to be of the same form with a 
lowered exponent, i.e., d k gi/dr] k oc r j L + n - l /' 2 - k (1 + o(l)). 



Functions with asymptotically polytropic behavior are generated by asymptotically 
polytropic distribution functions: 

Definition. A distribution 4>{£) (which is at least C°) is said to possess asymptotically 
polytropic behavior as £ — > 1, if there exists R 9 n > 1/2 such that 

</>(£) = <Ml-£) n - 3/2 + (62) 
as £ — * 1; here, 0_ > is a constant. 

Lemma A. 3. Distributions with asymptotically polytropic behavior generate functions 
gi{rj) with asymptotically polytropic behavior. 



Proof. By Corollary IA.21 g^ £ C^ +1 , since (f> is continuous. To prove the claim of 
the lemma we have to show that gi (and analogously its derivatives) behave according 
to ((Ml)- To that end set </>(£) = 0i(£) + </> 2 (£) with 



M£) = 0-(i - £) n " 3/2 , 2 (f) = - f ) n ~ 3/2 , (63) 
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where (f(£) = o(l) as £ — * 1. By definition, 



9i (v)= (M£) + M£)) H £2 - e ' 2v ) d£ = (M£) + M£))K l d£, (64) 



where we have set K = (£ 2 — e" 2r >)/2. The first integral can be computed explicitly by 
using hypergeometric functions; series expansion of the result leads to 



f (t)i(£)K t d£ = <p„ f (1 - £) n ^l 2 K l d£ = const rf+ n ~ l l 2 (1 + o(l)) 
as 7] — > 0. For rj small enough, the second integral can be estimated as 
C 4> 2 {£)K i d£ = f tf (£){l-£) n - 3 / 2 K i d£ 

< const (max{(p(£) \ £ € [e"", 1]}) rf^" 1 / 2 = o{rf+ n - 1/2 ) 



(65) 



Combining the results shows that 

gi (ri) = const 7] i+n - l/2 (1 + o(l)) . (66) 

By Corollary El 9i G C [i]+1 - To prove the lemma it remains to investigate the deriva- 
tives of gi. \ii> 0, the derivative of gi(rj) is given by (]51[). i.e., 

= i e - 2 " = const r? i+n - 3 / 2 (1 + o(l)) , (67) 

(XT? 

which shows that the derivative is asymptotically polytropic (with a lower exponent, 
as required). By iteration we obtain that all derivative behave correctly. If i = 0, we 
have 

which is again the required behavior. If % < 0, is differentiable only if </> satisfies 
the additional Holder continuity condition of Corollary IA.21 Also in this case it can 
be shown that the derivative behaves in an asymptotically polytropic manner; we omit 
the details here. □ 

Remark. If <j>(£) G C k (k > 0) and |(1 - £) k ^ k \£)\ < const (1 - £)~ 1+e , then higher 
derivatives of gi(rj) exist, cf. Corollary IA.2I In analogy to Lemma TA.3I it can be proved 
that these higher derivatives exhibit the correct asymptotically polytropic behavior 
provided that (1 — £) k (j)^ k \£) is asymptotically polytropic like <f>(£) in (|62[). 

Lemma A. 4 (The asymptotically polytropic case). Let <p(£) 6 C°(0, 1) be an 
asymptotically polytropic distribution function, see (|62 j) . (In the case I < —1/2, we 
require in addition that (f>(£) be Holder continuous with index a > —I — 1/2.) Then 
9i+3/2i 7 l) i s C 2 an d asymptotically polytropic as rj — > ; and 

lim o-(ry) = (69a) 
77— >o 
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Proof. Since (p G C°, Corollary IA.2I implies that ^+3/2(^7) is C 2 (where the Holder con- 
tinuity condition is needed in the case I < —1/2). Asymptotically polytropic behavior 
of (and its two derivatives) follows from Lemma lA. 31 Accordingly, we have 

log 91+3/2 = const + (I + n + 1) log 77 + log(l + o(l)) , (70) 

where the derivative of o(l) is of the order o(rj~ 1 ) as r/ — > 0; hence, from (|59[) . 

o-=[2l + 3 + (l + n + l)7 1 - 1 (l + o{l))]~ 1 ^0 (r/^0). (71) 

Furthermore, Eq. (|60ft results in 

^ = -a 2 [-(/ + n + 1)^(1 + (l))]^_L_ fo-0), (72) 

which proves the claim. □ 



B Global existence and uniqueness of solutions 



If a' (as a function of $7) is in C 1 (0, 1), the r.h.s. of the dynamical system (|23|) is C 1 ; 
existence and uniqueness of solutions follows directly from the theory of dynamical 
systems. These issues are less trivial when only continuity is imposed on a'. In the 
following we establish (global) existence and uniqueness of solutions of the dynamical 
system (|2*3*|) on the state space X = (U,V,Q) G [0, l] 2 x (0,1) under the assumption 
that a' be continuous on 

(0,1). 

We begin by introducing a new independent variable ( according to 



dX 



H' 



[(l-U + aU)VY 



(73) 



Then equation (|23cj) becomes d£l/dQ = —0(1 — Q)F(D,), where F(Q) is a positive 
function bounded away from zero for all O G [0,1]. Therefore, this equation has a 
unique global (in £) solution for any initial value $7(0) G (0, 1). In particular, f2(£) is a 
monotonic function mapping the interval (— 00, 00) to (0, 1). Eqs. (j23a() and (|23b|) take 
the form 



dU_ U(l - U) 
~dC ~ V(l -U + aU) 
dV 1-V 
~dC ~ 1-U + aU 



[3 + 21 - 2(2 + Z)Z7](1 - V) - 

(2U - 1)(1 — V + 2aV) + -^H 

dr\ 



1 + .-^) 77 
drj 



(74a) 
(74b) 



where a and a 1 are now regarded as functions of £. For any initial data (U(0), V(0)) G 
[0, 1] x (0, 1] the above system has a unique (local in £) solution. This follows immedi- 
ately from the Picard-Lindelof theorem, since the r.h.s. is continuous in £ and smooth 
in U and V. Translating back to the original variables, we have achieved existence and 
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uniqueness of solutions in the space [0, 1] x (0, 1] x (0, 1); a priori, existence of solution 
is local in A. 

Solutions with initial data in the space [0, 1] x (0, 1] x (0, 1) cannot reach the side face 

V = in finite A-time. This is a direct consequence of the fact that V~ 1 dV/d\ is 
bounded away from zero as V — > (where we use that U -f* 1/2). On the side face 

V = itself, existence and uniqueness of solutions is obvious, since the induced system 
has a simple structure. Combining these facts we have established (local in time A) 
existence and uniqueness of solutions on the state space [0, l] 2 x (0, 1). 

To obtain global existence of solutions we proceed as follows. Suppose that a (local) 
solution has been extended to a maximal interval (A_, A+) of existence. Assume that 
A+ oo. Since the solution ceases to exist for A > A+, it must leave the state 
space as A — > A+, i.e., O(A) — > as A — ► A + (note that O must be decreasing). 
However, O(A) satisfies the differential inequality dVt/dX > — FCl > — const O, therefore 
O(A) — ► const > as A — > A+, which is a contradiction. Hence, A+ = oo. Since we can 
show in an analogous way that A_ = — oo, the solution exists globally for all A. 



C Center manifold analysis 

Assume that a' (as a function of O) is in C 1 [0,1) so that we have an asymptotically 
polytropic C 1 -differentiable dynamical system on the state space [0, l] 2 x [0, 1). Then 
the following Lemma holds: 

Lemma C.l. The fixed point B\ = (0, j^p-, 0) is an attractor for interior orbits if and 
only ifa' < ^J. 

Proof. If erg < 1/(4 + 21), the fixed point B\ is hyperbolic: two eigenvalues are negative, 
one is positive. Analyzing the directions of the associated eigenvectors we find that there 
exists a one-parameter family of orbits that converges to B\ as A —* oo. 

Now assume that a' > j^l- In this case the fixed point B\ is a hyperbolic saddle (see 
Table ^) with the O = plane as unstable manifold; it immediately follows that B\ 
attracts no interior orbits. 

The case o~' = ^ s more delicate to analyze, since one of the eigenvalues is zero 
in this case; we make use of center manifold theory in our analysis. We define new 
adapted variables according to 

s:=(4 + 2J)u, y := {A + 2l)u + v- (4 + 2/) , z := rj ; (75) 

using these variables transforms the dynamical system 1)18(1 to a system of the form 

^ = N x (x,y,z) , ^ = y + N y (x,y,z) , ^ = -z + N z (x, y, z) , (76) 



where 

N x (x,y,z) = x 



3 + 21 - a' x - (1 + a(z) - a'(z)) (l + <r' a(z)x) [\ + y 



o 
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The center manifold can be represented by the local graph (x,h y (x),h z (x)), which 
satisfies the system of equations 

d x hy(x) < ^(x, h y (x), h z {x)) = ^(x, h y (x), h z (x)) , (77a) 

dx dz 
d x h z (x) — (x, h y (x),h z {x)) = — (x, h y (x), h z {x)) , (77b) 

and the tangency conditions hi(0) = 0, d x hi(0) = (i = y,z). By using expansions of 
the functions we are able to solve the system up to any desired order x n ; we obtain 
h z (x) = and h y {x) = —^^x 2 — ^^^^-x 3 + 0(x' 1 ). The center manifold reduction 
theorem states that the full nonlinear system is locally equivalent to the flow of the 
reduced system given by dx/d£ = N x (x, h y (x) , h z (x)) , dy/d£ = y, dz/d^ = —z, i.e., 

%-», Pa>) 

dz 

J!—. P8<=) 

The interior orbits are given by x > 0, y G K and z > 0; clearly, the flow of the reduced 
system prevents such orbits from converging to B\ so that B\ does not attract interior 
orbits. This concludes the proof of the lemma. □ 

The proof of Proposition 13.61 is now straightforward: 

Proposition EOfl Assume that there exists a variable u> such that c / | r? ( w ) is C 1 ^, e) 
(for some e > Q). If 

o~' > ^_^_2i ' or ' eau ^ va ^ en ^Vi n < 3 + 1 , (79) 
then the uj -limit of every interior orbit lies on (O = 0) D (V = 1). 

Proof. The case a' Q > 1/(4 + 21) (44> n < 3 + I) follows readily from Proposition 13.31 If 
cjg = 1/(4 + 21) (44> n = 3 + 1), there are a priori two attractors (which follows from 
analyzing the flow on = with the aid of the monotonic function Z, cf. H31|) ): a fixed 
point on (O = 0) fl (V = 1) or the fixed point B\. By Lemma IC. II the latter possibility 
is excluded, which proves the proposition. □ 



D The monotonicity principle 

In this section we prove an extension of the monotonicity principle presented in |16j ^ 
(see also [Til]). 

Consider an open set S C W 1 and its closure S. Let s be a point in S. We denote by 
B t (s) an open ball with radius e and center s, and by B e (s) its intersection with 5, i.e., 
B £ (s) = B e (s)DS. 
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Let Y : S — * (0, 1) denote a continuous function. For every s£5we define 

y_(s) := lim inf Y(x) , Y + (s) := lim sup Y(x) ; (80) 

the resulting interval [Yl(s), Y+(s)] we call the "limit set" of Y(x) as x — > s: 

'limy(x) := [Y^(s),Y + (s)] . (81) 

X— >s 

Note that Yl can be any number from to 1 (provided that Y is onto). 

Remark. Suppose that the set S possesses a "regular shape", i.e., assume that for each 
s E S there exists 5 > such that B e (s) is connected for all e < 5. (This is the case, 
e.g., for convex sets S). Under this assumption, for each a E 'lim-^g Y(x) there exists 
a sequence (x n ) ng N, x n £ S Vra, such that a = ]im n ^ too Y(x n ). This fact justifies the 
notation 'lint^s Y(x) for the interval [Y_(s), Y+(s)]. 

Remark. The set 'lim^^j y(x) depends in a continuous way on s: if ( 'lim x _»g y(x)) n 
[a — 5, a + 5] = for some a S K and <5 > 0, then for every 5' < 5 there exists a 
neighborhood of s such that ( 'linx^^ Y(x)) n{a — 5' ,a + 5'] = for all a G U (s). 

Remark. For a function Y that possesses the range (—00,00) the limit set 'liml" is 
defined in complete analogy; in this case, 'limy is a subset of [—00,00]. Note that Y 
can be manipulated by using a diffeomorphism between (0, 1) and (—00,00) to obtain 
the original set-up. 

Consider a dynamical system 

§ - '<«> < 82 > 

that is defined on a future invariant open set S C M n ; in the case when / is merely 
continuous, we require existence and uniqueness of solutions of (|%2"|) to hold. We say 
that f(x) does not diverge as x approaches the boundary OS of S, if the condition 
±00 'lim^—jg f( x ) holds for all s € dS. (In fact, in the following theorem we will 
meet the slightly less restrictive assumption that ±00 'lim^j f(x), Vs G dS\(dS)- m f.) 

Suppose that there exists a C 1 function Z : S — > K that is strictly monotonically 
increasing along all orbits in S. We define the subset (35)j n f of dS as 



(dSkf = <s £ dS I 'lim Z(x) = {inf Z} } . (83) 

[ x-^s S J 

Based on the definitions and concepts introduced above we are now able to state and 
prove the following theorem: 

Theorem D.l (Monotonicty principle). Consider a dynamical system dx/dt = 
f(x) on the future-invariant open set S C M n ; f(x) is required to be non-divergent 
on dS\(dS)- m {. Suppose that there exists a C 1 function Z : S —> M that is strictly 
monotonically increasing along orbits in S, whereby the derivative Z = X7Z ■ f(x) is a 
function Z : S — > (0, 00). Then for all x £ S, 

uj(x) C J s G dS\(dS) mi I ( 'lim Z(x) 3 supZ) V ( 'lim Z{x) 3 0) 1 , (84) 

I X — >S g x — *s \ 
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Proof. Since Z is strictly increasing along all orbits in S, the infimum inf Z (and the 
supremum supZ) is not assumed in S. For each x G S we have Z(x) > inf Z, and 
hence a; (a;) (<9«S)i n f- 

Assume that there exists x G 5 such that 

w(z) n j s G dS\(dS) ioi I ( 'lim Z(x) ^ sup Z) A ( 'lim Z(x) $ 0) 1^0, (85) 

; 1 2!^S £ X^S 7 J 

and let s be an element of this set. Consequently, there exists a sequence (t n ), n G N, 
such that i n — ► oo as ti — > oo and limn^oo (x) = s (where <^ denotes the flow of 
the dynamical system). Since Z is monotonically increasing along orbits, the sequence 
Z((f>t n (x)), n G N, is strictly monotonically increasing. The limit of this sequence 
is finite, lin^-nx, Z((f) tn (x)) < oo, because 'lining Z(x) ^ sup s Z (and 'lim x ^ s Z(x) 
is a closed interval). In contrast, the limit hin^^oo Z((f) tn (x)) need not exist; since 
'lim Z{x) j$ 0, there exists e > such that Z(<pt n (x)) > 2e for sufficiently large n G N. 

a;— >s 

Similarly, by assumption, there exists F G (0, oo) such that \f((f>t n (x))\ < for 
sufficiently large n. 

Choose a neighborhood U(s) of s such that 

'lim > e and 'lim |/(^ n (x))| < F (86) 

for all a G £7(s)- Let d be defined as 2d = in£ s eS\u(s) |s — s\, i.e., 2d is the minimal 
distance between s and S\U(s). For n large enough, \<fit n (x) — s\ < d; hence the time 
r it takes the flow to transport <t>t n {x) out of U(s) can be estimated by r > d/F. (We 
show below that the flow must indeed transport <pt n {x) to S\U{s) again, at least for 
infinitely many n.) Since Z > e during the time r we obtain 

Zifa+M) > Z(<t>t n+T {x)) > Z(c(> tn (x))+ed/F . (87) 

for an infinite number of values of n. This is a contradiction to the finiteness of 
lim^oo Z((f) tn (x)). 

It remains to show that the flow cj) t transports 4>t n (x) out of a (sufficiently small) 
neighborhood U (s) of s at least for infinitely many n. Assume, indirectly, that the 
statement were false. Then, for all neighborhoods V(s) there exists some m G N such 
that 4>t m +t{x) G V(s) for all t > 0, which means that (j>t(x) — ► s (t — > oo). Consequently, 
since lim^oo Z((f) t (x)) is finite, the sequence Z{<p t {x)) must contain a subsequence that 
converges to zero. This is a contradiction to our assumption 'lining Z(x) ^ 0. □ 
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